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ABSTRACT 

We present the results of a series of observations of the near- and mid- infrared polarisa- 
tion properties of the Seyfert 2 galaxy NGC1068. Our data agree well with previously 
published results in showing the need for a separate polarisation mechanism in the 
near infrared apart from scattering. We find that the predictions of a simple model in 
which this component arises through absorptive dichroism due to aligned grains within 
the extended warm (~ 400K) dust fits the data reasonably if the obscured background 
source is itself due to dust emission (at temperature > 1000K). By considering the 
change of polarisation with wavelength we show that the extinction to this hot dust 
region is in the range Ay = 20 — 40. Consideration of the observed data then leads us 
to the conclusion that if viewed face-on, NGC1068 would have a strong near-infrared 
excess similar to Seyfert 1 galaxies. Comparison with other independent measures of 
the extinction to the active nucleus itself lead us to the conclusion that the hot dust 
must provide screening equivalent to at least Ay = 40, and possibly much higher. We 
speculate that this component alone may be the 'classical' torus discussed in terms 
of the unified model, and the more extensive mid-infrared emission may arise from 
circumnuclear molecular cloud material, and dust in the ionisation cones. 

Key words: galaxies: individual (NGC1068) - galaxies: Seyfert - galaxies: nuclei of 
- infrared: galaxies - galaxies: active - polarisation 



1 INTRODUCTION 

The unified model for Seyfert galaxies proposes that all 
types of Seyfert galaxy are basically the same, but that 
a dusty molecular 'torus' obscures the broad-line region 
(BLR) in many systems. The classification into Seyfert 1 
or 2 (Seyfert 1: broad permitted lines and Seyfert 2: narrow 
permitted lines) then depends on the inclination angle and 
hence on whether we can see the BLR. Although the con- 
cept of the unified model existed before the discovery by An- 
tonucci & Miller (1985) of broad Ha emission in spectropo- 
larimetry of the nearby luminous Seyfert 2 galaxy NGC1068, 
that result was seen by many as the major confirmation of 
the likely validity of the model. Although other Seyfert 2's 
have also been found to have broad permitted lines in scat- 
tered light, this is by no means universal. Indeed Heisler, 
Lumsden & Bailey (1997) found a strong correlation be- 
tween the 'warmth' of the IRAS 60/im to 25/^m ratio and the 



ability to detect such scattered broad-line radiation. Further 
studies that enhance our understanding of why some Seyfert 
2s appear consistent with the unified model and others do 
not are clearly therefore of value. 

In this context, infrared polarimetry of AGN provides 
valuable extra constraints on the polarisation mechanisms 
active in Seyfert 2's. Since the unified model infers the pres- 
ence of a dusty torus obscuring the AGN core, we might 
expect that longer wavelength polarimetry will be able to 
probe deeper into regions within the plane of the torus. 
Hence infrared polarimetry may let us see scattering in 
Seyfert 2's which is otherwise shielded from view. In addi- 
tion, many Seyfert 2's with generally 'warm' IRAS colours 
are known to have polarisation rising into the near infrared 
(Young et al. 1996a). The mechanism generally invoked to 
explain such behaviour is linear absorptive dichroism, in 
which aligned dust grains preferentially absorb one plane of 
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polarisation. When such dust screens a bright background 
source, the transmitted light is therefore polarised. As shown 
by Hildebrand (1988), observations of this aligned dust com- 
ponent at wavelengths where it emits as well as at wave- 
lengths where it absorbs can be used to study the nature 
and geometry of the dust that so heavily obscures the AGN 
along our line of sight. 

Since it is already known from the spectropolarimet- 
ric observations that NGC1068 contains a hidden type 1 
core, has relatively high polarisation for a Seyfert 2 (Heisler, 
Lumsden & Bailey 1997), and has rising infrared polarisa- 
tion (Young et al. 1995), it makes a perfect target in which 
to test these ideas. Its proximity (14.4 Mpc; Tully 1988) and 
luminosity also make it an ideal candidate. 

There has already been considerable previous work on 
infrared polarisation in NGC1068. The most extensive early 
broad-band polarimetry data are reported by Bailey et al. 
(1988), who found a strong rising red continuum in the near- 
infrared polarised flux and evidence for a gradual change in 
the position angle of polarisation between the optical and 
near infrared. They explained these as a combination of po- 
larisation by scattering and due to transmission of a strongly 
reddened background source through aligned dust grains. 
They also found a much more marked change in polarisation 
angle (~ 70°) between 4 and 5/im due to the same aligned 
grains being seen in emission at longer wavelengths. This 
large swing in polarisation is seen again in the spectropo- 
larimetry of Young et al. (1995) who also included polari- 
sation by aligned grains in their model, having found that 
scattering fitted the data well below 1/xm but was deficient 
beyond that wavelength in explaining both the polarised flux 
and the position-angle data. 

In addition, Aitken et al. (1984) presented Wfim spec- 
tropolarimetry, showing that the percentage polarisation re- 
mained constant through the 9.7/im silicate absorption fea- 
ture. This showed that the cause of the lO^im polarisation 
could not be absorption by aligned dust grains, since that 
would have resulted in a rise in polarisation through the 
absorption feature. However, they were unable to decide be- 
tween the other alternatives: emission from aligned, optically 
thick (since the silicate absorption feature seen is relatively 
weak) dust grains or polarised emission from the active core 
itself. The position angle of polarisation seen at 10/um is 
consistent with that found by Bailey et al. (1988) at 5/im. 

A consistent picture has begun to emerge from all these 
results. The optical data near the nucleus can be explained 
by a combination of a dominant electron-scattering com- 
ponent and a lesser dust-scattering one. The near-infrared 
polarisation is a combination of electron scattering and the 
dichroic component modelled by Young et al. (1995). The 
longer wavelength data is then consistent with emission from 
aligned grains, since the emission and absorption polarisa- 
tion position angles from the same grains are roughly or- 
thogonal. 

However, two aspects were clearly lacking in all of this 
work. Since much of the activity in NGC1068 occurs within 
1-2 arcseconds of the nucleus, high spatial resolution ob- 
servations are required to separate the different polarisa- 
tion mechanisms that may be present. Some high-resolution 
imaging polarimetry has been published by Young et al. 
(1996b) but this was obtained in moderate seeing condi- 
tions. The only other imaging data that exists in the near 



infrared is that of Packham et al. (1997), but the resolution 
of this data is > 1.2 arcseconds. All other previously pub- 
lished data were acquired using aperture bolometer devices. 
The second major lack is imaging polarimetry in the mid- 
infrared. Without such data it is not possible to relate the 
observed near-infrared data to the behaviour at longer wave- 
lengths with a high degree of confidence. We have therefore 
gathered high-resolution (~ 0.5 arcsecond), complete near- 
and mid-infrared imaging polarimetry of the central regions 
of NGC1068 using the AAT. In section 2 we give details of 
these observations, in section 3 we present the results of our 
observations, in section 4 we analyse these results in terms 
of possible polarisation mechanisms, in section 5 we discuss 
the dust distribution around the nucleus of NGC1068 and 
finally in section 6 we present our conclusions. 



2 OBSERVATIONS 

All of the data presented in this paper were obtained on the 
nights of 9, 10 and 11 August 1995 and 11 October 1997 
at the Anglo Australian Telescope. The 10/im data were 
acquired with the ADFA mid-infrared camera, NIMPOL 
(Smith, Aitken and Moore, 1994), which at the time used a 
128 2 SiGa array, with a pixel scale of 0.25"/pfxel. A cold wire 
grid and a warm rotating CdS half-wave plate were used as 
the analyser for the polarimetry. We used a standard broad 
band 8-13/im filter for the observations, and chopped and 
nodded ~ 20" to sky so that the object was always present 
on the array. Extended emission at 10//m is sufficiently weak 
(eg Telesco & Decher 1988) that this poses no problems for 
the nuclear polarimetry. Approximately 7 hours of actual 
on-target data were acquired. The image is diffraction lim- 
ited, and all data taken during non-photometric conditions 
have been discarded. We observed NGC1068 at airmasses 
ranging between 1.2 and 1.8. 

The near-infrared data were obtained on the nights of 
11 August 1995 and 11 October 1997 with the common-user 
camera IRIS, which uses a 128 2 HgCdTe array. A warm ro- 
tating half-wave plate and a cold beam-splitting Wollaston 
prism were used to obtain the linear polarimetry: the Wol- 
laston allows for very efficient polarimetric observations, and 
also allows reliable data to be taken even in non-photometric 
conditions. Again the pixel scale was 0.25"/pixel. The seeing 
was estimated at 0.6" during these observations. After each 
cycle of four waveplate positions we nodded the telescope to 
a separate sky position since the galaxy is bright in the near 
infrared beyond the nucleus. Approximately 20 minutes of 
data were obtained on source at each of J, H and K n . Con- 
ditions were photometric for the H and J band observations, 
and partly photometric for the K n observation. The photom- 
etry given in Section 2 is derived solely from those frames 
taken when conditions were photometric. Unfortunately, the 
conditions meant that it was difficult to obtain good 'sky' 
frames to subtract from the data at K n due to clouds reflect- 
ing thermal emission from the ground. The result is a higher 
than expected background noise level in the polarisation im- 
ages. We therefore also observed NGC1068 at H and K n on 
the night of the 11th October 1997 with the same setup to 
check our results. We obtained 16 minutes of on-source data 
at K„ and 8 minutes at H on this night. The seeing was 
good for the K n data (~0.75 arcseconds), but poorer for the 
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H band data (~1 arcsecond). Conditions were photometric 
during this night. This is evident in a considerable decrease 
in the noise in the polarised flux images. The polarimetry 
derived from both dates is consistent within the errors. 

Flux calibration for the mid-infrared data was obtained 
from matching observations of the mid-infrared standards 
BS6832 (r? Sgr; S(10jum) = 197 Jy), BS8636 {(3 Gru; S(10^m) 
= 933 Jy) and BS7525 (S(10pm) = 93 Jy), and has an es- 
timated uncertainty of 10%. Flux calibration for the near- 
infrared 1995 data was derived from SA241-251 and for the 
1997 data from HD24849 (Carter & Meadows 1995). 

The polarimeters used with both IRIS and NIMPOL are 
very stable, and instrumental effects result in errors as low 
as 0.1% in polarisation and less than 1° in position angle. In 
order to check the performance of the system with IRIS we 
also observed standard stars of known polarisation, as well 
unpolarised standards and stars observed with a wire grid 
(of known polarisation properties) in the beam. This allows 
the zero-point of the position angle to be set. The resulting 
errors found were in line with the expected performance. 
For NIMPOL, the polarisation standard was the BN object 
in Orion, assumed to have a polarisation angle of 118° at 
10/xm. 

The near-infrared data were reduced as follows. The 
data were flat-fielded using dome flats. Offset sky frames 
(several arcminutes from the source) were obtained using the 
same setup as for the object data. The frames for each half- 
wave plate position were median-filtered to remove back- 
ground sources. The resultant four sky frames (one for each 
half-wave plate position) were scaled to the median level of 
the sky within each group of four individual object frames, 
and the result subtracted from the object frames (this tech- 
nique proved resilient to residual sky variations from tests 
carried out by subtracting sky from individual sky frames 
themselves, which could otherwise provide an unwanted ad- 
ditional background 'noise' signal in the final polarisation). 
The resultant images were then registered and combined into 
separate mosaics for each half-wave plate position. These fi- 
nal mosaics were combined to form Q and U Stokes images, 
and hence polarisation maps. This procedure is valid be- 
cause of the dual-beam nature of the instrument. Variations 
in sky brightness, or transparency, are correctly allowed for 
because the frames are scaled when combined to form the 
Stokes parameters ensuring that the total flux is conserved 
in any image. 

The raw 10/im data is stored with the chop already 
subtracted but the corresponding nod position is removed 
later in software. The result is then properly sky subtracted. 
Flat fields were made from pairs of sky observations, one at 
high (~ 1.6) and one at low (~ 1.0) airmass, which were 
subtracted and normalised. The sky-subtracted, flat-fielded 
data were registered to remove any small shifts in the object 
position with time, since the telescope was tracking without 
guiding. The shift observed is only ~ 1 arcsecond per hour. 
Frames in which the background increased by more than 
50% compared to others in a given sequence were discarded. 
Such frames were generally affected by cloud. The final data 
were then grouped into Stokes Q and U parameters. Note 
that since NIMPOL uses a wire grid analyser and not a 
Wollaston prism, residual variations in atmospheric trans- 
parency can affect the final results. We tested for this by 
breaking the data into sequential blocks and compared the 



results. We found no evidence that cloud affected any of the 
data used in the final images. 



3 POLARISATION MAPS AND 
PHOTOMETRY 

3.1 The near-infrared data 

The near-infrared polarisation maps are shown in Figure 
1 and the corresponding polarised- flux images (being flux 
multiplied by degree of polarisation) are shown in Figure 2. 
We have derived an error estimate for the polarised flux from 
the scatter in the counts well away from the nucleus, and 
in Fig. 1 only plot polarisation points which lie more than 
3cr above the mean background level. The images shown are 
derived from the 1997 data for H and K n , and from the 1995 
data for J. It is worth noting that the centroids of the direct 
flux and polarised flux are identical at all three wavebands 
to within less than one tenth of a pixel (0.02 arcseconds). 

The measured polarisation is given in Table 1. The er- 
rors are derived by measuring the polarisation on every indi- 
vidual subset of four half-wave plate positions. We note that 
the measurement in the small aperture is highly sensitive to 
factors such as seeing, telescope drift etc, and therefore has 
a much higher error than the other values quoted. The re- 
sults are in reasonable agreement with both Bailey et al. 
(1988) and Packham et al. (1997). We believe that the dif- 
ferences found by Packham et al. between their data and 
that of Bailey et al. may be largely due to factors such as 
seeing variations. The core of the emission from NGC1068 
is clearly strongly polarised, but the profile of this core is 
narrow, so the resultant measured polarisation in a small 
aperture is strongly dependent on factors such as seeing. 

The most immediately obvious features present in Fig- 
ure 1 are the centro-symmetric vector patterns characteristic 
of scattering seen both northeast and southwest of the core. 
This provides a direct constraint on the geometry of the 
central regions of NGC1068. Since two scattering cones are 
seen in the present data, but only one in in the optical and 
UV, the counter cone must be hidden by an extended screen 
at shorter wavelengths. The only likely candidate for this 
screen is the disc of the galaxy. This implies that the cones 
in NGC1068 are definitely not aligned along the plane of the 
galaxy, in agreement with previous evidence that the north- 
eastern cone illuminates the near side of the disc whereas 
the southwestern cone illuminates the far side (cf Figures 1 
and 2 of Bland-Hawthorn et al. 1997). In particular, maps 
of the [OIII] emission show that the southwestern emission 
region vanishes behind the larger-scale molecular ring that 
surrounds the nucleus at a distance of about lkpc, whereas 
the northeastern emission is clearly visible. The same phe- 
nomenon can also explain the asymmetry in the X-ray emis- 
sion (Wilson et al. 1992). 

The polarisation vectors are clearly evident spread in 
an arc between the common axis of the large-scale ionisa- 
tion cones and radio emission (eg. Wilson & Ulvestad 1987, 
Evans et al. 1991; position angle ~ 30°) and that of the inner 
radio jet (eg. Ulvestad et al. 1987; position angle ~ 10 — 30°). 
Whether or not this is one single continuous structure is not 
absolutely clear from the present data, however, which hint 
at there being two preferred directions in which the scatter- 
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ers lie, coincident with the two axes of the large- and small- 
scale radio emission. Data from HST may be invaluable in 
determining the true situation. 

Our data are in good agreement with the images pre- 
sented by Young et al. (1996b). The 'shadow' they noted in 
their H band data is also present in our data. We plot the 
image profile along the 45° direction (ie along the long axis 
of the images in Fig. 1) in Figure 3 for J, H and K. Note 
that our sign convention is opposite to that of Young et al., 
since positive offsets here are in the same sense as in Figure 
1. In Fig. 3 we have used the H band data from 1995, when 
the seeing was better than for the 1997 H band data. The 
K band data from both 1995 and 1997 are consistent, but 
we show the former to ensure consistency with the other 
plots. The clear dip in the profile seen at J and H at an 
offset of ~ 1 arcsecond agrees well with the H band profile 
plotted by Young et al. along a slightly different position 
angle. They attributed this to the possible shadowing of the 
scattering region by the torus. This is consistent with the 
presumed orientation of the ionisation cones if the torus is 
also aligned in the same sense. 

Photometry in the three wavebands from both the Au- 
gust 1995 and October 1997 data are given in Table 2. Errors 
are estimates based on the scatter between the frames ob- 
tained. The larger error on the K n data from August 1995 
reflects the fact that not all of the data were photometric. 
For the 1997 data, we used all of the K n data but only the 
first two minutes of the H band data since the seeing deterio- 
rated after that point (clearly visible from the change in the 
individual image profiles). We have used only data taken in 
seeing better than 1 arcsecond. Because of the limited field 
of view, it is difficult to judge if the sky has been fully cor- 
rected for in these images. Light from NGC1068 dominates 
the images at all wavebands. We have adopted an average 
correction for the sky based on the detected counts 10 arc- 
seconds from the nucleus. If this is an overcorrection, the 
tabulated fluxes will increase by less than 5%. 

We can compare these results with those of Glass 
(1995), who published photometry over a twenty- year times- 
pan in a 12 arcsecond aperture. There are several prelimi- 
nary corrections that need to be made, however. Glass used 
a standard single-element photometer so that only 12 arc- 
second aperture data is published. Because of the mask used 
in the imaging polarimetry our data is restricted to a 6 arc- 
second aperture at most. Since NGC1068 varies slowly with 
time in the near infrared (Glass 1995), we cannot compare 
our data directly with other existing photometry. However, 
we can use such photometry to derive an estimate of the flux 
within a smaller aperture from the Glass data. We compared 
the Glass data at Julian date 2446314 with that from Bailey 
et al. (1988) in a 6-arcsecond aperture at JD~ 2446450. This 
is a sufficiently small separation in time that the variability 
should be negligible. We ascribed the difference to emission 
in the annular ring between 3 and 6 arcseconds radius from 
the peak. As noted by Glass, the source that varies has the 
colours of hot dust: it is essentially the same as the bright 
core seen within 2 arcseconds of the nucleus in our images. 
It is likely, therefore, that the annular ring has little or no 
variability in flux, since it lies outside the hot-dust emission 
region. Assuming this is correct, we can subtract the derived 
annular flux from later data from Glass (JD 2449615) which 
is closer in time to our August 1995 data. Using this approx- 



imation we derive a 6 arcsecond aperture flux density from 
the Glass data of 761mJy for K, 363mJy for H and 198 mjy 
for J. The errors on these numbers are entirely limited by 
the systematics of the approximations made, but it is likely 
that these values are reliable to within 10%. The results are 
in good agreement with our own 1995 data. 

We can also compare our photometry with that of Pack- 
ham et al. (1997). Their H band result would agree with ours 
if we made no additional sky correction. Their K band result 
is discrepant with our data even if we make this assumption 
(their result being ~ 400mJy brighter than the value given 
in Table 2). 

Lastly, we note that the 1995 and 1997 data are in rea- 
sonable agreement with each other, and provide no strong 
evidence for a large change in brightness at these wave- 
bands over this time period. Glass (1997) has reported 
that NGC1068 may now be declining in brightness at near- 
infrared wavelengths, although the magnitude of the de- 
crease is small and roughly consistent with our own results. 



3.2 The mid-infrared data 

The lO^im image is shown in Figure 4. The effective diffrac- 
tion limit of the AAT at 10/im is ~ 0.6", and the dust emis- 
sion is clearly not a point source. The data also clearly show 
that the polarisation structure arises from an extended re- 
gion (all points with a signal-to-noise ratio greater than 5 
in the polarisation are shown). In addition, there is a trend 
for the polarisation to increase away from the peak of the 
flux. The position angle of the polarisation is essentially con- 
stant across the region within the errors at 55° ± 10°. The 
measured aperture polarisation is given in Table 1. 

Aitken et al. (1984) obtained aperture spectropolarime- 
try at 10/im of the nucleus of NGC1068, showing a relatively 
featureless polarised flux spectrum. They concluded that the 
likeliest cause was emission from warm aligned dust grains, 
but could not rule out emission from the active nucleus it- 
self. However, their data conclusively rule out an origin for 
the polarisation in absorption by dust since this would give 
a peak in the percentage polarisation through the silicate 
absorption feature at 9.7/im. The result they derived was 
1.39 ± 0.09% polarisation at a rather poorly defined posi- 
tion angle (measured on two separate occasions at 44.1±3.2° 
and 59.4±2.2°in a 4.2 and 5.6 arcsecond diameter aper- 
ture respectively). Our observed position angle is consistent 
with the latter, though the measured polarisation is slightly 
higher in the 4.5 arcsecond aperture. 

Dust emission can be seen up to 5" from the nucleus, 
with the extended emission largely aligned with the ob- 
served ionisation cones (the inner high contours are also 
aligned with this direction as noted previously by Braatz 
et al. 1993). The measured flux density within a central 6 
arcsecond aperture is 24.7±2.1Jy, and within a 10 arcsecond 
aperture 26.0 ± 2.2Jy. The surface brightness at the lowest 
contour in Fig. 4 is ~ lOOmJy/sq. arcsecond. 
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4 ANALYSIS 

4.1 Polarisation mechanisms and the location of 
the AGN core 

Young et al. (1995) found the dominant contribution to the 
polarisation at short wavelengths (A < 1.5/im) is from scat- 
tering (mostly from electrons). We can check this by ex- 
amining the deviation from the expected centro-symmetric 
scattering pattern. As shown by Capetti et al. (1995), if 
only scattering is present then the average of the differences 
between the observed position angle in any pixel and the 
expected position angle (which is just the normal to the 
vector joining that pixel and the central source from which 
light is scattered into our line of sight), should be zero. Of 
course, this is a semi-iterative process since the exact loca- 
tion of the scattering centre is a priori unknown. However, 
it is trivial to iterate towards a solution for both the source 
of the scattered light and the deviation from the expected 
centro-symmetric pattern. 

Figure 5 shows the observed departure from a pure- 
scattering vector pattern for all three wavebands. It is clear 
that the deviation in the expected vector pattern reaches 
a maximum at position angle ~ 115° (ie for vectors with 
predicted position angle of 30°) at all wavelengths. The ob- 
served pattern is consistent with that expected if polarised 
flux with a constant position angle were superposed on the 
observed scattered flux (see Section 4.2). Indeed this pattern 
is clear from Figure 1 as well. Only by masking out a region 2 
arcseconds in radius around the central nucleus do we actu- 
ally see no deviations from the expected scattering pattern. 
It is clear therefore that another mechanism contributes to 
the observed polarisation at all wavebands. 

We masked out all the data within 2 arcseconds of the 
peak when deriving the actual scattering centre because of 
this additional mechanism that contributes near the nucleus, 
as well as those data with signal-to-noise ratio in the po- 
larised flux less than 3<r above the mean background. We 
found we could not determine the location of the scattering 
centre reliably by projecting orthogonal vectors back from 
the polarisation vectors to find a common centre. Residual 
background results in a 'blurring' of the polarisation posi- 
tion angles, so that some of the points of intersection lie well 
away from the nucleus. We therefore calcuated, for each pixel 
in turn, the difference between the observed position angle 
and the position angle expected if that pixel were the cen- 
tre. We then found the global minimum of this function, and 
this position is the scattering centre. This is far more stable 
to the effects noted, but at the cost of giving a less precise 
solution than would be possible with the other method given 
better data. This technique gave offsets from the respective 
flux centroids of 0.11 ±0.18, 0.09 ±0.21 and 0.34 ±0.39 arc- 
seconds at J, H and K„ respectively, which are all formally 
consistent with no offset. 

We note that Marco, Alloin & Beuzit (1997) found the 
centroid of the K band flux tentatively aligned with the lo- 
cation of the UV scattering centre (Capetti et al. 1995), and 
the radio sources SI and S2 (Gallimore et al. 1996), from as- 
trometry obtained using adaptive optics and offsets between 
the peaks of the optical and IR light. Our results are there- 
fore entirely consistent with their analysis. We note that 
their result depends on determining a difference between the 
position of the peak of the I band light and the peak of the 



K band light. Although it might seem unusual that the peak 
of the J band light should therefore be coincident with the 
peak of the K band light (since I is closer to J in wavelength 
than J is to K), the reason is simple. The J band centroid 
is dominated by the same polarised core that is seen at H 
and K even though the contrast with the electron scattering 
cones is much less. The same feature is unlikely to appear 
at I, as is clear from Section 4.2. It is the centroid of this 
feature that determines the location of the AGN core as we 
show below. 

Lastly, we consider the factors that may affect our so- 
lution. These are systematic global position angle errors, 
random position angle errors and the role of seeing. Sys- 
tematic global errors in the position angle calibration are 
constrained to be < 1° from observations of polarised stan- 
dard stars, which results in an entirely negligible change in 
the location of the scattering centre. We checked for effect of 
random position angle errors using a Monte-Carlo analysis 
of our observed data, drawing a new set of position angles 
from the set by randomly shifting the angle consistent with 
the observed error in the position angle. This had a negligi- 
ble effect on the location of the centre. The effect of seeing is 
to blur the scattering pattern, especially near the bright nu- 
clear continuum source. We have attempted to quantify this 
by simulating a simple ionisation cone irradiated by a cen- 
trally peaked Gaussian source, with noise levels throughout 
consistent with our data. This pattern was then smoothed 
to mimic the observed seeing, and rederived the polarisation 
from the model Stokes Q and U images. The result shows 
that seeing effects the observed pattern only in the inner 
0.5 arcseconds around the flux centroid. Since we exclude 
this data anyway, seeing clearly has little affect on our re- 
sult either. Therefore we can be confident that our result is 
correct, and that the scattering centre is coincident with the 
K band flux centroid. 

4.2 The nature of the polarised core 

As indicated by the analysis above, there is at least one other 
polarisation mechanism present within ~ 1 — 2 arcseconds 
of the peak of the flux. Young et al. (1995) originally pro- 
posed that this was absorptive dichroism on the basis that 
no scattering law appeared to fit their spectropolarimetry. 
That suggestion was consistent with the results of Bailey et 
al. (1988) who also concluded that the cause of the ~ 70° 
change in position angle between the near infrared and the 
mid infrared is the switch from emission to absorption by 
aligned grains as wavelength decreases. A lower limit to the 
temperature of the dust providing the screen can be set by 
the observation of Bailey et al. (1988) that the position an- 
gle of polarisation changes somewhere between 4 and 5/im, 
from which we conclude that T dust > 350K. 

The fact that the change in polarisation position angle 
is not closer to 90° may indicate one of two possibilities: first 
that the background source seen at A < 4/u,m contributes to 
some extent at 10/im too but is polarised itself with a dif- 
ferent position angle (so at short wavelengths, the position 
angle is determined by the combined effects of dichroic ab- 
sorption and polarised emission from the background, and 
at 10/im from emission from both); second, that there is yet 
cooler dust that is still shielding the 10/im emission, so that 
the observed position angle is a blend of the emission and 
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absorption processes (one test of this may be obtained by 
determining the polarisation in the sub- mm regime) . 

All our data are consistent with the background source 
being hot dust with temperatures near the sublimation point 
(see below), where grains are unlikely to be aligned by a 
magnetic field. Therefore if the grains are aligned by a mag- 
netic field, the first possibility is ruled out. Alternative mech- 
anisms such as alignment due to radiation pressue in the 
ionisation cone may give rise to a similar pattern (eg radi- 
ation pressure: Dopita et al. f 998) but it is difficult to see 
how these can give rise to alignments in both hot (> f 000K) 
and warm/cool (< 400K) dust components. It is also known 
that NGC1068 has a prominent silicate absorption feature 
at 9.7/im (Roche et al. 1984), indicating the last explanation 
may be the most probable. If it is true, we would predict that 
the observed position angle may change through the silicate 
absorption feature from the 45° measured at 5/im by Bailey 
et al. to the ~ 55° we measure. Unfortunately the data of 
Aitken et al. (1984) are not sufficiently accurate to measure 
such small shifts. 

We have modelled the wavelength dependence of the 
observed polarised flux in the core. The model is simple 
and assumes only three parameters that are allowed to vary. 
These are the temperature of the dust, Tdust, the emissivity 
of the dust (assumed to vary as a power law, A - ™), and the 
extinction to this dust source as measured in the K band, 
Ak- The polarised flux is assumed to be represented by an 
optically thin grey body with this temperature and emis- 
sivity, scaled by an appropriate Scrkowski polarisation law 
(Serkowski, Mathewson & Ford 1975 and references therein) 
and then reddened using a standard A -1 ' 75 extinction law. 
The Scrkowski law has the form 

p(X)/Pmax = exp [-iHn 2 (Xmax/^j] 

where we have adopted 'typical' values for moderately ex- 
tinguished sources in our galaxy (Whittet et al. 1992) of 
K — 1.15 and \ ma x = 0.6/im (the effect of changing X ma x 
and K is small for the wavelength range of interest). 

The results are indicative only given the few data points 
we are trying to fit. It may be possible, given small aperture 
spectropolarimetry between 1 and 4/xm, to provide better 
discrimination between the models. However, a typical 'good 
fit' to the data can be obtained with the following param- 
eters: T dust = 1200l2oo, n = 1-5 ± 0.5 and A K = 1.5 ± 0.5 
(or Ay = 17 ± 6). The errors given here are not formal 
uncertainties but represent the range of values that pro- 
duce reasonable fits to the three data points. The results 
are also correlated in the sense that larger n gives lower val- 
ues of Tdust, and larger values of Ak require larger values 
of Tdust. However, it is clear that an extinguished hot dust 
source can explain the observed data. Since we might expect 
the hottest dust present to be near the sublimation temper- 
ature, the 'likeliest' parameters are those with high Tdust 
and A K (ie T du st=1500K and A K = 2). 

We also considered a simple power law for the under- 
lying source rather than a grey body. A reasonable fit can 
again be derived with similar Ak assuming the underlying 
source varies as A, which, however, is inconsistent with typ- 
ical synchrotron spectra. For the dust temperature derived, 
there is essentially no flux below 1/im, explaining why and 
J and K n centroids are similar to each other, and not to I 
band data. 



We can estimate what fraction of the total flux at J 
within a 3 arcsecond aperture arises from the dust assuming 
the model with both large T dust and Ak is correct. If we use 
the result of Thatte et al. (1997) that dust contributes 90% 
of the direct flux at K, and our model above, we predict that 
dust will contribute ~ 50% of the total flux at H, and ~ 10% 
at J. This is in agreement with our findings in Section 4.1 
that even at J there is a source of polarisation other than 
scattering, and that the J flux centroid is coincident with the 
K flux centroid, since the centroids at both wavebands are 
determined largely by the centrally peak dust component. 
The same hot dust contributes similar flux density at 2 and 
10/im, or ~ 5% of the total at 10/im. Since in reality there 
will be a continuous range in dust temperature with distance 
from the source, it is likely that up to 25% of the 10/im 
omission can come from dust hotter than 1000K. 

Lastly, we can obtain an estimate of the position angle 
and dimensions of the inner regions of the obscuring ma- 
terial from our data. The position angle of the symmetry 
axis of this 'torus' is ~ 30° (90° from the deviation in the 
scattering pattern shown in Figure 5) . This agrees well with 
previous estimates of the position angle of the torus (Young 
et al. 1995, 1996b, Miller et al. 1991), and of the large-scale 
radio emission (Wilson & Ulvestad 1987). The observed data 
are also marginally suggestive of an extended source aligned 
perpendicular to this, since the profiles along this direction 
and along the cone differ in extent by 0.1 ± 0.05 arcseconds. 

To derive better estimates of the size of the resolved 
inner-torus region we created a simple model in which a 
polarised core was superposed onto background scattered 
emission. The relative strengths of the core and background 
were taken from the actual K„ data. The polarisation was 
also taken from the observed data. We then convolved this 
model with the observed seeing. The resultant polarisation 
maps, though highly idealised, could be compared with the 
data to determine the actual extent of the polarised core. 
The results indicate that this region has approximate ex- 
tent 7 x 15pc. We show in the next section that this ex- 
tent is consistent with the observed mid-infrared emission 
from the ionisation cones. Clearly, high spatial resolution 
(~ 0.1 arcseconds or better) IR imaging polarimetry that 
can sample the fine structure in this region is desirable to 
estimate the true size of the torus. 



4.3 Mid-infrared emission 

As shown by both Braatz et al. (1993) and Cameron et al. 
(1993), it is possible to relate the observed dust temperature 
and its distance from the heating source to the source lumi- 
nosity, after making suitable assumptions about the grain 
properties and assuming that the dust sees the source di- 
rectly). Cameron et al. derive the following equation, as- 
suming an emissivity proportional to A~ 1,5 : 

( R \ = ( 0.05/im x L \° B / T dust y 2 - 7 
\0.14pcy \a x 1.5 x 10 11 L© / U500K7 

Here a is the grain size in microns, L is (effectively) 
the luminosity of the core and Tdust the dust temperature. 
Clearly, for fixed R and L, the observed dust temperature is 
weakly inversely proportional to the grain size. Therefore, we 
might expect that it is the physically smallest grains (with 
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radii~ 0.005/im) that give rise to the hottest dust emission 
at any point and we adopt this factor in the following anal- 
ysis. 

First we consider the dust in the ionisation cone. Al- 
though we do not have an estimate of the dust temperature 
in the ionisation cones, we can assume a lower limit, which 
then gives a lower limit to the luminosity. The fact that 
we see dust emission at all at 10pm implies T dust > 150K, 
just from a consideration of the sharpness of the drop in 
flux density with wavelength of a suitable grey body. If this 
dust is radiatively excited by the nucleus, this implies that 
the dust in the ionisation cones sees a luminosity of at least 
~ 3 x 10 1 Lq . This is a factor of two higher than the esti- 
mated bolometric luminosity along our line of sight and may 
indicate anisotropic emission. This is not a new conclusion; 
Miller et al. (1991) found that the properties of the elec- 
tron scattering region were also consistent with anisotropic 
emission from the core, as did Young et al. (1997). 

Returning to the issue of the compact hot-dust emis- 
sion discussed above, we see that for the same parameters, 
dust can be heated to 600K at a distance from the core of 
8pc. Note we assume a lower limit to the temperature for 
the marginally 'extended' emission at K, since the 1500K 
dust visible at all three near infrared wavebands is likely 
to be an unresolved point source in our data (there is no 
visible extent at H and J so we cannot apply this analy- 
sis). Given the many assumptions that have gone into this 
simple analysis, the fact that this distance agrees with the 
limits set by our data is highly encouraging. The more ex- 
tended dust around the nucleus should be dealt with using 
a proper radiative transfer model, since the source is likely 
to be screened by the hotter dust in that case (cf Pier & 
Krolik 1993, Efstathiou et al. 1995). It is interesting to note 
that the same parameters predict a dust temperature of ap- 
proximately 100K at the radius of the large scale molecular 
ring seen in CO band images of NGC1068, which may just 
be warm enough to explain the excess 10pm emission seen 
in that ring coincident with the ionisation cones (Bland- 
Hawthorn et al., 1997). 

We modelled the observed 10-pm polarisation structure 
as a combination of a point source (a reasonable approxima- 
tion for the hot dust) and an extended component (warm 
dust). We assumed that the point source was completely un- 
polarised. The same results are obtained if the point source 
is polarised along a different axis to the extended component 
however. We convolved the point source with a suitable Airy 
function and subtracted this profile from the 10-pm data un- 
til the percentage polarisation was constant across the nu- 
cleus. The best results are obtained if we allow the point 
source to have peak height 25% of the observed peak flux. 
The resulting polarisation pattern after subtraction then has 
p~ 1.4%. This result is also consistent with that expected 
from the model outlined in Section 4.2, where we found that 
between 5 and 25% of the emission at 10pm could arise in 
dust at temperature greater than 1000K (ie dust where the 
alignment is likely to be destroyed). 



5 COMPARISON OF MID AND NEAR 
INFRARED DATA 

Young et al. (1995) modelled the polarised core of NGC1068 
as arising due to absorptive dichroism, with an extinction to 
the emitting source of Av ~ 45. However, this was based on 
the assumption that the contribution of the dust to the con- 
tinuum emission seen in their 3" x 3" aperture is ~ 15% 
at H, whereas Origlia et al. (1993) estimate this contribu- 
tion is ~ 30% (in a slightly larger aperture) and Thatte 
et al. (1997) derive an estimate of 65 ± 5% within a similar 
sized aperture. At K, Thatte et al. find that dust contributes 
89 ± 5% of the light within the same aperture, which is also 
larger than the value given by Origlia et al. (70%) within a 
slightly larger aperture. It is worth comparing these values 
with the crude estimates derived in Section 4.2. There we 
found ~ 50% of the total H band light should be due to dust 
which is midway between the Thatte et al. and Origlia et al. 
estimates, but given the crudity of our procedure probably 
consistent with both. Packham et al. (1997) also estimated 
the required Ay to explain the observed data using a similar 
method to that given below, though placing the emphasis on 
the residual polarised flux after scattering had been allowed 
for, and derived Av ~ 35. It is useful to consider this prob- 
lem independently of the 'best fit' scattering model to see 
what limits can be placed. 

The shape of the near-infrared spectrum is very sensi- 
tive to the dust temperature, assumed emissivity and extinc- 
tion to the emission region as shown in Section 4.2. As noted 
there, changing any of these parameters has an effect on the 
others since they are correlated in any fit. It is therefore 
very helpful to derive independent estimates of any of these 
quantities. One method of deriving an independent estimate 
of the extinction is to compare the relative polarisation at 
10pm and K. The dust that shields the hot emission source 
evident in the K band is clearly cooler (and more spatially 
extensive from our 10pm image) as noted above. As shown 
by Hildebrand (1988), the intrinsic percentage polarisation 
seen due to aligned dust grains both in emission and absorp- 
tion can be related according to 

Pern = Pabs/T\(abs), 

if the result is taken at the same wavelength for both. This is 
clearly impractical, and a useful relation between absorptive 
polarisation at K and emissive polarisation at 10pm must 
depend on the wavelength dependence of the dust cross- 
sections. 

There are several other key assumptions in deriving a 
useful result from this simple relationship. First, it is as- 
sumed that scattering is not important at short wavelengths. 
This is true if we use a sufficiently small aperture near the 
peak of the flux at K„, and at H to a lesser extent. It also as- 
sumes that the 10pm emission is optically thin which is prob- 
ably a reasonable approximation (note the fact that there is 
an observed silicate absorption feature as shown by Roche 
et al. (1984) does not violate this condition). The same 
grains must contribute at both wavelengths, and the wave- 
length dependence of the absorption/emission cross-sections 
are taken from the observed dependence of the known ex- 
tinction law, from which tk ~ 1-7tlo, and hence the same 
for the cross-sections. Therefore, we apply a correction of 
0.6 to the observed K n band polarisation. 
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We cannot use the observed values directly in this for- 
mula, since it relates solely to the properties of the dust. We 
can, however, correct these values, by estimating what frac- 
tion of the polarised light is due to scattering rather than 
dichroism, and what fraction of the direct light is due to stars 
rather than dust emission. For the former, we have derived 
aperture photometry on the polarised flux images within 
an aperture of diameter 3 arcseconds. We subtract off the 
extended contribution underlying the central source, assum- 
ing it is due to scattering (since this is where the scattered 
broad permitted lines arise). The result is assumed to be the 
polarised flux due to dichroic absorption of a background 
source. Of course, if the 'background correction' is not rep- 
resentative of the scattering properties near the core (for 
example, if the electron scatterers themselves are sharply 
peaked near the core) this will be an overestimate of the 
polarised flux due to dichroism. We can be confident this is 
not the case, however, just from the observed behaviour of 
the position angle of polarisation near the core. This essen- 
tially shows that within the innermost 1 arcsecond, scatter- 
ing cannot dominate the light. The derived corrections are 
small for the region within 1.5 arcseconds of the core, and 
we find that 90 — 95% of the polarised flux in this region 
arises in the strongly peaked feature seen in Figure 3. 

For the stellar contribution to the direct light, we use 
the values given by Thatte et al. (1997) but note that there 
appears to be an inconsistency between the results of Thatte 
et al. and Origlia et al. (1993). The difference in the result 
is small for the K-band data (Ay is ~ 3 magnitudes larger) , 
but significant if we use the H-band data. As a result we 
do not consider the H-band data in detail, but note that 
a similar result (in terms of Av) to that derived at K is 
obtained if we use the lower Origlia et al. value. 

We also need to correct for possible extra sources that 
contribute to the 10pm polarisation. As noted above it is 
possible that the hot dust contributes 25% of the light but 
none of the polarisation at 10/im, and we have removed a 
point-source contribution at this level. The main uncertainty 
in the mid-infrared polarisation, however, is the dust opacity. 
Assuming the emitting regions are partially opaque implies 
that some of the emissive polarisation is itself absorbed as 
noted before. Hence the observed polarisation is actually 
reduced since the net effect of the absorption is to reduce 
the contrast between the axes of maximum and minimum 
polarisation (which is also the likely reason why the position 
angles are not orthogonal as noted in Section 4.2). Therefore, 
the result derived here should be thought of as an upper 
limit. 

Therefore, assuming that p a t s > 1.4% and using p em ~ 
4.7%, we derive tk < 3.3, or Ay < 37 assuming a standard 
IR extinction law. This sets the upper limit to the values the 
data allow. Since the estimated extinction to the AGN itself 
is > 80 (Jackson et al. 1993), the cool dust is clearly only 
providing a fraction of the observed extinction. Our result 
is much larger than the observed optical depth in the sili- 
cate absorption feature (Roche et al. 1984 derive Ay ~ 8), 
but reasonably consistent with the C-H band stretch ab- 
sorption seen at 3.4/im (Bridger et al. 1994 find Av ~ 22) 
the fits derived for the near-infrared polarised flux in Sec- 
tion 4.2 (where we found Av ~ 23 a reasonable fit) and the 
value derived by Packham et al. (1997) (A v ~ 35). The dis- 
crepancy with the silicate absorption feature is most likely 



explained by infilling from the hot dust emission washing 
out the absorption feature, as shown by full radiative trans- 
fer models (Efstathiou et al. 1995, Pier and Krolik 1993). 
Clearly therefore, if the hot dust is cospatial with the AGN 
core, then that dust itself must provide the bulk of the ex- 
tinction to the core along our line of sight. 



6 CONCLUSIONS 

We have obtained new near- and mid-infrared imaging po- 
larimetry of NGC1068. Our data are consistent with other 
observed results, but have generally higher spatial resolu- 
tion than anything yet published. We can see clearly from 
our data that scattering alone fails to fit the observed data 
at all wavebands, and not just at K as noted previously by 
Young et al. (1995). We have shown how this deviation can 
be explained by absorptive dichroism of a background hot 
dust source, and find that at least part of the hot dust must 
have a temperature > 1000K. The screening dust must have 
temperature > 350K, since the observed change in position 
angle expected when moving from aligned grains emitting 
to aligned grains absorbing arises at ~ 4 — 5pm. It is likely 
that even cooler dust is also present which is shown only by 
the fact that the change in position angle in the infrared is 
not exactly 90°. 

We have derived independent estimates of the extinc- 
tion to the background hot dust source seen through the 
obscuring screen. The results, from fitting the observed near- 
infrared polarisation data and from comparing the near- 
and mid-infrared polarisation give results that are consis- 
tent with Av — 20 — 40. This is considerably less than the 
estimate for the extinction through the observed molecular 
material to the AGN core itself (eg Jackson et al. 1993 find 
that Ay = 80 from consideration of the optical depth ob- 
served in the circumnuclear molecular gas) . The large visual 
extinction to the hot dust however explains why previous 
authors have noted the absence of a near-infrared excess 
in NGC1068 as typically seen in Seyfert l's (eg Edelson & 
Malkan 1986, Cameron et al. 1993). If the hot dust were 
unobscured, the combined spectrum of it and the more ex- 
tensive warm dust would indeed peak near 5pm. It may be 
that the hot dust component is actually the 'classical' torus, 
whereas the warm dust is distributed more widely and dif- 
fusely throughout the circumnuclear molecular clouds and 
the narrow line region. This is in line with the results of 
Heisler et al. (1997) who found a strong correlation between 
the size of the electron scattering region and the obscuring 
material in a sample of far-infrared selected Seyfert 2's. 

It is clear from our results that near-infrared polarime- 
try offers a potentially valuable insight into the nature of 
the obscuring material. In particular, it allows us to deter- 
mine the properties of the warm dust that is nearest the 
AGN core. The clearest need for the future is high spatial 
resolution spectropolarimetry to define the polarisation be- 
haviour as a function of wavelength, and for polarisation 
predictions to be included in radiative transfer modelling of 
the dust absorption/emission near AGN. The former will be 
especially useful when applied to techniques such as that dis- 
cussed in Section 4.2, as well as greater information of the 
overall polarisation properties of dust absorption features, 
which may show weak features in the change in position an- 
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gle not previously considered. There is also clear scope for 
these methods to be applied to other nearby Seyfert 2s with 
warm far-infrared colours. 
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J band polarisation 
Aperture Diameter Polarisation Position Angle 



(arcsec) (percent) (degrees) 

2 2.25 ±0.25 106.8 ±2.0 

4.5 1.59 ±0.15 105.3 ±1.0 

6 1.35 ±0.10 103.9 ±1.0 
H band polarisation 

Aperture Diameter Polarisation Position Angle 

(arcsec) (percent) (degrees) 

2 3.87 ±0.20 117.0 ±2.0 

4.5 3.06 ±0.10 115.0 ±1.0 

6 2.58 ±0.05 114.5 ±0.5 
K n band polarisation 

Aperture Diameter Polarisation Position Angle 

(arcsec) (percent) (degrees) 

2 4.57 ±0.5 120.2 ±2.0 

4.5 4.19 ± 0.08 120.3 ±0.5 

6 4.07 ±0.08 119.6 ±0.5 
N band polarisation 

Aperture Diameter Polarisation Position Angle 

(arcsec) (percent) (degrees) 

2 1.30 ± 0.05 49 ±3 
4.5 1.80 ± 0.05 57 ±6 

Table 1: The measured polarisation within a circular aper- 
ture around the peak of the observed total flux. 
Aperture Diameter Flux Density 
(arcsec) (mJy) 
J band photometry 

3 93±10 
6 198±10 

Aug. 1995 H band photometry 

3 243±10 

6 390±15 
Oct. 1997 H band photometry 

3 240±15 

6 420±15 
Aug. 1995 K n band photometry 

3 650±30 

6 800±30 
Oct. 1997 K„ band photometry 

3 553±20 

6 785±10 



Table 2: The measured photometry for NGC1068 in a cir- 
cular aperture around the peak of the flux. 
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Figure 1: Near infrared polarisation maps of the nucleus of NGC1068. The K n data are shown at the top, H in the middle 
and J at the bottom. A vector 1 arcsecond long represents 5% polarisation. The contours are of the total flux and are scaled 
arbitrarily. 
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Figure 2: Near infrared polarisation images of NGC1068. The K n data are shown at the top, H in the middle and J at the 
bottom. The scale is the same as in Figure 1. 
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Figure 3: Cross section through the images shown in Figure 2. The K n data are shown at the top, H in the middle and J 
at the bottom. The dip in the profile at —1 arcsecond offset at H and J is the same as the feature reported by Young et al. 
(1996) 
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Figure 4: Mid infrared polarisation map of the nucleus of NGC1068. The data in the greyscale image have been smoothed 
slightly to enhance the contrast in the faint extended emission. The data are displayed on a power-law scale (surface brightness 
to the third power). The contours are evenly distributed in on the same scale, but are otherwise arbitrary. 
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Figure 5: Deviation of the observed position angle of polarisation from that expected from a centro-symmetric scattering 
pattern. Pure scattering would be represented by a straight line through zero. Again, the K„ data are shown at the top, H in 
the middle and J at the bottom. The greatest departure from zero is always at approximately 115°. The width of this feature 
in the J band data is a reflection of the relative weakness of the additional polarisation component in this case. 
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